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Abstract 



The necessary and sufficient condition of separabihty of a mixed state of 
any systems is presented, which is practical in judging the separability of a 
mixed state. This paper also presents a method of finding the disentangled 
decomposition of a separable mixed state. 
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Entanglenment is an important element both in quantum theory [|I| and in quantum in- 
formation 0-^. A mixed state p, which acts on Hilbert space Hi®H2, has no entanglement 
or is separable iff p can be written as the following form 

i 

where p\ and pf are states of Hilbert space Hi and respectively, and > 0, X^Pj = 1. 
A separable system satisfies the Bell inequality, but the converse is not necessarily true [0 . 
Peres has shown in Ref. that a necessary condition (NC) for separability of a state p is 
the positivity of its partial transformation p^ . Here is defined as 

(2) 

It has been shown that for 2x2 systems (the dimensions of Hi and H2 are 2 and 2, 
respectively) and 2x3 systems, the condition is also a sufficient one 0. A necessary 
and sufficient condition(NSC) for m x n systems with the aid of the concept of a map 
has been presented by Horodecki et.al in Ref. 0, but it is difficult to judge whether a 
state is separable or not in practice with this method. W.Diir et.al [l^ show a NSC of 
a certain family of mixed states. In this paper, we shall present a NSC of separability 
of p for any system. Our main idea is based on the idea proposed by Wootters in Ref. 
p!T|JT^ . Compared with Horodecki et.al s work, the NSC presented in this paper is more 
explicitly algorithmic, so it is practical in operation. The organization of this paper is as 
follows. First, we collect some facts related to our result. Second, the NSC of any mixed 
state is presented. Meanwhile, how to get the disentangled decomposition of a separable 
state is also presented. Finally, we present an example of a bound entangled state as an 
illustration. From this example, one can know how to judge whether a mixed state is 
separable or not. 

In order to get the results, we collect the facts as follows: 

l.Any mixed state p of m x n systems which has / nonzero eigenvalues can be de- 
composed into k pure states. And the number of possible decompositions of p is infinite. 
These decompositions can be obtained by some transformation Ukxi whose columns are 
orthonormal vectors {k is greater than or equal to /) [|ri[]. For example, 

I k 
i=l i=l 

where \xi) , unnormalized, is a complete set of orthogonal eigenvectors corresponding to 
the nonzero eigenvalues of p, and Xj) is equal to the ith nonzero eigenvalue. Then the 
decomposition \zi) of p can be given by 

I 

= '^Uij\xj) , i = l,2,---,k (4) 



where the vectors \zi) are not necessarily orthogonaL 
2. A pure state of an m x n systems is denoted as 

1^) = All 111) + Ai2 |12) + ■■■ + Ain\ln) + ■■■ + Ami + ■■■ + A^n \mn) 



(5) 



where |11) , ■ ■ ■ |1?t,) , ■ ■ ■ |ml) , ■ ■ ■ \mn) are the standard basis vectors of the m x n system. 
We introduce the "vector" Aj 



Ai — (v4ji, Ai2, ■ ■ ■ Ai 



1,2, 



m 



(6) 



It is not difficuh to find that the pure state is separable if and only if all the "vectors'' 
Ai are parallel with another , which is denoted as 



Ai Ma II ■■■ II Am. 
Note that the relation(0) is true if and only if the following is true 



(7) 



{Aii,Ai2) II (Ai, A22) 

(Aii,Ai3) II (^21,^23) 



{Ami, Am2) 
{Ami: Am'i) 



^ {All, Am) II (A2I, A2n) II ■ ■ ■ II {Ami: Amn) 

That is to say there exist (m — l)(n — 1) parallel pairs. For example, having parallel pairs 
{Aii,Ai2) II (A2i,A22),(Aii, A2) II (Ai,A2) means that Ai/Ai = A2M22, AiMsi = 
^i2M32, respectively. 

3. We define ^^'") and with respect to each parallel pair as follows: 



B' 1^*) , Cr = {^\ ^'■) , r = 1, 2, ■ ■ ■ , (m - l)(n - 1) 



(9) 



where is a mn x mn symmetric transformation matrix in the standard basis and there 
are only 4 nonzero elements of B^. For example, for (^11,^12) || (^21,^22), [-S^]^„_|_2 ~ 
[i?"'^]^^2 1 = ~'^'[B^]2n+i ~ \-B^]n+i2 ~ other elements of [B^] are all zero. It is 

obvious that |ci| = if and only if {Aii,Ai2) \\ {A2i,A22)- Similarly, each parallel pair in 
the relation(H) can be expressed by its |cr| being equal to zero and this is also the reason we 
choose the elements of the matrix B'' . For 4x4 systems, there exists nine transformation 
matrices as follows: 



{AIUA12) 


II (Ai, A2) - 


-IB'] 


1,6 = 


-1J5']2,5 = 


1, [B\2 = 


1. [B%,i = -1, 


{Aii,Ai2) 


II (^31,^32) - 




1,10 ~ 


-^AB%,,= 


- 1> [B%,2 = 


= lJ5']io,i = -1, 


{Aii,Ai2) 


II (Ai, A2) - 




1,14 ~ 


-^AB\is 


= l>[5']l3,2 


= 1,[B\,^, = -1, 


{Aii,Ais) 


II (^21,^23) - 




1,7 = 


-1, [^^3,5 = 


1^ [B%,s = 


1, [B%,i = -1, 


{Aii,Ais) 


II (^1,^3) - 




1,11 = 




- 1> [B%,s = 


= lJ5']ii,i = -1, 
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(^1,^3) 
(^1,^14) 
(^11,^14) 
(^11,^4) 



(^1,^3) 
(^21,^24) 
(^31,^34) 
(^41,^44) 



J3,13 



1,8 ~ 

1,12 = [B 



-1, [B' 



lJ^']l3,3 
4,9 = 1 J^\4 = 1 J5%2,1 



-1, 
-1, 
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1, [B'' 



1, 



1,16 J4,13 J13,4 l"^ Jl6,l 

the other elements of [B^'] (r = 1, ■ ■ ■ , 9) are all zero. Similarly, for higher dimension 
systems, we choose B^ of each parallel pairs so that each Cr is zero. 

4. From the above facts, one can get that any state p is separable iff there is a pure 
state decomposition oi p, p = J2i l^i) {^il with each \zi) being separable, i.e., each pair in 
Eq.(|D is parallel for all \zi) . 

In the following, we will discuss the NSC of separability of m x systems. Firstly, we 
will prove the following lemma. 

Lemma. For any density matrix p of mxn systems in Eq.(^, the statement {An, A12) \\ 
{A21, A22) is true for each pure state \zi) in Eq.(H) if and only if 



(10) 



1=2 



where /' < / {I, I' is the number of eigenvalues of p, r^(r^)* respectively.) and AjS, in 
decreasing order, are the square roots of eigenvalues of the matrix r^(r^)*, where t^j 



{xi \ x^) and \xi) is defined in Eq.(|D and 



x: 



m 



Eq.®. 



Proof: Obviously is a symmetric matrix. For the symmetric matrix 
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one can 



find a / X / unitary matrix [uij] which can diagonalize the 



in the following way 



1 T 

VT V 



(11) 



where \yi) = J2j ^ij , v^j is the complex conjugation of Vij, and AjS are the square 
roots of eigenvalues of matrix r^(r^)* [|ri|]. Though is dependent on the choice of 
eigenvectors \xi) of p, Aj is not. Every decomposition oi p, p = X^i {zi \ , can be obtained 
by ki) = EjKi \yi)- Here 



0216 



16*2,1 



0-2,26 



401 



102,2 



ai,«e' 

^2,^6 



afc,ie 



ak,2e 



(12) 



whose columns are orthonormal vectors and there is no loss of generality in taking 
each ttij to be real and nonnegative. If /' < I, one can add (Z — /') dummy states 
\yi'+i) ) \yi'+2) , \yi) being equal to zero vectors. Thus it can be obtained that \zi) = 
YjjUij \xj), where = p*!/* [there may exist many transformation matrices which can re- 
alize this kind of decomposition. The set which contains all these transformation matrices 
is denoted by f/^xi)- If (^11,^12) || (^21,^22) for each \zi) , from fact 3, one has 
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zl) = 0, (13) 

where {zj) = \z*):, then one has 

<A = (E 4 cos 2%A,)2 + iJ2 4 sin 2%A,)2 

i=2 j=2 

<(E<A-r- (14) 
i=2 

Because of the positivity of Aj and aj^, one can get 

I' 

J=2 

and then, Ai = Eti a^Ai < ELi E5'=2 = Si=2 ^i' 

I' 

a' = Xi-J2Xi<0. (16) 
Here we have used Eq.(|15|) and Y.i=i o-a = 1- 

Conversely, if Ai < Yli=2 -^j? one can find a kind of decomposition of p, p = J2i \ 
[see Appendix], where {zi\ zj) = 0, so that (An, A12) || (^21,^22)- Therefore the lemma is 
proved. 

Similarly, from the proof of the Lemma one knows that for an m x n system the 
NSC for each of the (m — l){n — 1) pairs being parallel is < 0(r = 1, 2, ■ ■ ■ , (m — 
l)(ri — 1)). If each < 0, there exist (m — l){n — 1) sets of transformation matrices, 
^fcx/> ^fcx/> ■ ■ ■ > ^fc™^^'*^"~^^ 6^cl^ of which yields decompositions of p with the existence of 
the corresponding parallel pair. Thus we get the main result in this paper: 

Theorem. A NSC of separability of states of an m x n system is that the corresponding 

< (r = 1, ■ ■ ■ , (m — l)(n — 1)) and the intersection Ukxi of the (m — l)(n — 1) sets of 
transformations is not empty. 

Proof: The sufficient condition is obvious. The task left is to prove the necessary 
condition. If p is separable , there must exist a transformation matrix u which can transform 
the decomposition of p, p = Y.i=i l^i) {^il into p = Y.i l^i) {zi\ with the existence of (n — 
l)(m — 1) parallel pairs in each pure state \zi). From the Lemma, one can get each < 
in each \zi), and the matrix u belongs to the intersection Ukxi', i-e. Ukxi is not empty. Thus 
the theorem is proved. 

According to the above theorem, one may take the following steps (i)-(iv) to judge 
whether a state p of an m x n quantum system is separable or not. 

i). Calculate the nonzero eigenvalues and eigenvectors of p (see Eq.(H)). 
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ii) . Calculate (m— 1) matrices r^(r^)* according to B^{r = 1, 2, ■ ■ ■ , (m— 1)) 
and obtain the square roots of the eigenvalues of t'^{t^)* (see Eq.(|Tl])). 

iii) . Judge whether a^{r = 1,2, ■ ■ - , (m — l){n — 1)) is greater than zero or not (see 
Eq.(|TOD). If all a*" < 0, then go to step (iv), otherwise p is inseparable. 

iv) . Judge whether the intersection of {m ~ l){n — 1) sets of transformation Ukxi is 
empty or not. If the intersection is empty, p is inseparable, otherwise p is separable. 

The steps from (i) to (iii) are easy to follow, but step (iv) is very difficult to operate at 
first sight. It is natural to ask whether step (iv) can be done or not and whether the method 
is practical . To answer this question and illustrate the method, we give an example of a 
bound entangled state. 

A state p of a 2 X 4 system is 



P = o 



10 
10 
10 
1 
1 
10 
10 
10 



1 
1 




1 
1 

1 



(17) 



p has five nonzero eigenvalues tj(i = 1, 2, ■ ■ ■ , 5), 

_ 1 

_ _ _ 1 

t3 — — — -, 

The corresponding unnormalized eigenvectors \xi) can be 

= (0,0, 0,1/^,0,0, 0,0), 1x2) = (0,0,0,0,1/^,0,0,0), 



(18) 



'i, 0, 0, 0, 0, J^, 0, 0), 1x4) = (0, 0, 0, 0, 0, J^,0), 



ka) = 

ks) =(0,0, 

The three matrices B'^{r = 1, 2, 3) corresponding to the three parallel pairs are 



(19) 



w 8.0,0,0,0, 



(^11,^12) II (^21,^22) = B^ -^ = -1,-82,5 = -^5,2 = 1' 

{An, An) II (^21,^23) ^ Bl, = Bl, = -l,Bl, 



(^11,^14) II (^21,^24) -Bis = -^8,1 = -Ij-^is 



^5,3 
^5,4 



1, 
1. 



(20) 
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Then the symmetric matrices tI 
are 



{Xi\ x'j) = {xi\ B"" 



l,2,---,5,r = l,2,3) 



"o 





o" 




"o 








o" 





















|o 













1 

8 




|o 














1 

4 





,r^ = 








1 

8 
















1 

8 


o| 













1 

8 

















































1 

8 









(21) 



It is easy to get that = 0, < 0, < 0, so we must turn to step (iv). According to 
the theorem, there are surely three sets of transformation matrices Ukx5 each of which can 
result in the decomposition \zi) of p, p = J2i=i \zi) {zi\ with the corresponding pair being 
parallel . The m^xs can be written as 



MfexS 



an 



ai2 

0'k2 



ai5 

025 



(22) 



Each column in m^xs is an orthonormal vector. If three sets of Ukxb have intersec- 
tion, the state p is separable. This condition is equivalent to the existence of a set 
On, ■ ■ ■ , ai5, ■ ■ ■ , afci,...,afc5 satisfying the following equations 



\zi) = an \xi) + a,i2 \X2) H 

{z,\B^\z*) = Q, 
{z,\B^\zl) = Q, 
[{z,\B^\zl) = Q. 

From Eqs. (|21|) and (p3D, one can get that 



1,2, 



k) 



(23) 



*2 



(24) 



The complex number a*j can be expressed as a*j = hijC^^'^ where hij is a real and positive 
number. From the first equation of Eq.(^ 
Because Yh=i bjj = 1, the equality bi2 = bi^ 



one can get that J2i=i bi2bi4 = bf 



hi = hi 



•i2 



b. 



•i3 



■- bi4 can be obtained. Then Eg . (|24|) results in 

bi4 = = bi 



8ii+3ei2 a _ 0ii+9i2 n 



bi5 



-'11 
39ii+0i2 



(25) 



Because each column in Eq.(^2]) is orthonormal, one can obtain that 
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Eil 62e^(^«i-e>2) =0, 
< Eti 62e^(^.i-^«) = fofg.^i^^ ^ 0, (26) 

It is easy to find that Eq.(|2|) has no solution since Eti hfe'^^^^-^^''^ = 0, Eti 6- e^ ''"''^'""' = 
0, and Ei=i ^ = do not hold true simultaneously.. Therefore the intersection Ukx5 
of three sets of transformation matrix is empty and p in Eq.([l7D is inseparable. 

From the above example, one can see that for any mixed state the problem of finding 
the intersection of our sets of transformation matrices comes down to the problem of 
finding the solution of a set of quadratic equations with the condition of the solution being 
orthonormal vectors. The condition of orthonormality that makes the quadratic equations 
more easily solved. Especially, if there are not many nonzero elements of the symmetry 
matrix r^', it is very easy to solve the set of quadratic equations. In addition, the form 
of these equations is dependent on the choice of eigenvectors of p, which implies that one 
may simplify the form of these equations with appropriate eigenvectors of p. 

In summary, we have presented a new necessary and sufficient condition for separability 
of state of any system. According to the above theorem, (m — 1) x (n — l) conditions should 
be tested in order to judge whether a state is separable or not. Though in mathematical 
form this method is more complex than the one provided in [^], it is practical in judging the 
separability of a state. Meanwhile, our theorem provides a method to find the disentangled 
decomposition of a separable state. It is also worth mentioning that for 2 x 2 systems, the 
condition that the partial transposition(PT) is positive and the condition that each < 
in our method are equivalent. This suggests that there may be an essential connection 
between the PT of p and a^. 
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APPENDIX 

Considering an mxn system, a state p of this system, which has / nonzero eigenvalues, 
has a decomposition, p = J2i IVi) {Uil ; which satisfies {yi\ = Xi6ij, where = \y*) 
and is a transformation matrix corresponding to a parallel pair (r = 1, 2, ■ ■ ■ , (m — 
l){n - 1)). If Ai < Ei=2^i is the number of nonzero eigenvalues of r'"(r'")*), one can 
find a decomposition of p, 

1 ' 

l^^) = T7fE«ue*'M2/j)' ^ = 1,2,---,4A;, A; = l,2,3,--- (27) 
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where Ak is greater than or equal to /, \aij\ = 1 and their signs are decided by Fig. 1, 
\yi'+i) ) \yi'+2) , \yi) being zero vectors. Phase factors can be found to make J2j e^*^^ Aj = 



||Tl|, which guarantees {zi\ zl) = 0. Obviously, the matrices [aij] is a Ak x I matrix. 
From Fig. 1, one can find that 



04x4 



-1 +1 +1 +1 
-1 +1 -1 -1 
-1 -1 +1 -1 
-1 -1 -1 +1 



^8x8 



+1 +1 


+1 +1 +1 +1 +1 


+1 


+1 +1 


+1 +1 -1 -1 -1 


-1 


+1 +1 


-1 -1 +1 +1 -1 


-1 


+1 +1 


-1 -1 -1 -1 +1 


+1 


+1 -1 


+1 -1 +1 -1 +1 


-1 


+1 -1 


+1 -1 -1 +1 -1 


+1 


+1 -1 


-1 +1 +1 -1 -1 


+1 


+1 -1 


-1 +1 -1 +1 +1 


-1 



(28) 



According to the symmetry, one can find the higher-dimension matrix a^kxAk easily. If the 
number / of nonzero eigenvalues of p is less than 4k, one may take the first / columns of 
«4fcx4A:- For example. 



3x5 



+1 +1 


+1 


+1 +1 


+1 +1 


+1 


+1 -1 


+1 +1 


-1 


-1 +1 


+1 +1 


-1 


-1 -1 


+1 -1 


+1 


-1 +1 


+1 -1 


+1 


-1 -1 


+1 -1 


-1 


+1 +1 


+1 -1 


-1 


+1 -1 
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Figure Captions 

Fig. 1. The diagram showing the sign of the oefficient Ukj in Eq. (16). 
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FIGURES 




Fig 1 
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